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Abstract 

In this paper we calculate the turbulent heating rates in the so- 
lar wind using the Kolmogorov-like MHD turbulence phenomenology 
with Kolmogorov's constants calculated by Verma and Bhattacharjee 
[1995b, c]. We find that the turbulent heating can not account for 
the total heating of the nonAlfvenic streams in the solar wind. We 
show that dissipation due to thermal conduction is also a potential 
heating source. Regarding the Alfvenic streams, the predicted tur- 
bulent heating rates using the constants of Verma and Bhattacharjee 
[1995c] are higher than the observed heating rates; the predicted dis- 
sipation rates are probably overestimates because Alfvenic streams 
have not reached steady-state. We also compare the predicted tur- 
bulent heating rates in the solar corona with the observations; the 
Kolmogorov-like phenomenology predicts dissipation rates compara- 
ble to the observed heating rates in the corona [Hollweg, 1984], but 
Dobrowoly et al.'s generalized Kraichnan model yields heating rates 
much less than that required. 
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1 Introduction 



The proton temperature T of the solar wind, regarded here as a single mag- 
netofluid parameter, is observed to decrease slower that what adiabatic cool- 
ing would predict [Schwenn, 1983; Marsch et al, 1983; Schwartz and Marsch, 
1983; Gazis, 1984; Freeman and Lopez, 1985; Lopez and Freeman, 1986; Free- 
man et al, 1992]. These observations indicate that the protons in the solar 
wind is heated in its transit. Note that in this paper we will not discuss the 
temperature evolution of electrons and alpha particles. Also, since the solar 
wind has complicated flow behaviour, the simplistic arguments presented in 
this brief report are of a qualitative nature. 

Various attempts have been made to explain the observed heating. The 
proposed sources of heating are shocks [Whang et al, 1990], turbulence [Tu et 
al, 1984; Tu, 1987; Tu, 1988; Verma et al, 1995a], heat conduction [Gazis, 
1984], interactions with neutral particles [Isenberg et al, 1985] etc. Whang et 
al [1990] performed observational and simulational studies on the corotating 
shocks in the solar wind between 1 and 15 AU, and showed that in this region, 
shocks are major heating sources in the solar wind, with most of the heating 
confined to 1-5 AU. They found that the increase in the entropy per shock 
is approximately 0.8 x 10~^^ J/K/proton. Note that entropy is deflned as 
ks ln(T^-^/n), where T is the temperature, n is the particle density, and ks is 
the Boltzmann constant. They also derived that the temperature variation 
due to the observed entropy increase yields T{r) oc r~°'^^, where r is the 
distance from the sun. Note that without any heating, the temperature of 
the solar wind would be proportional to r~^/^. 

In a recent paper, henceforth referred to as paper I, Verma et al [1995a] 
calculated turbulent heating in the solar wind using the existing MHD tur- 
bulence phenomenologies, Kolmogorov-like MHD turbulence phenomenology 
[Marsch, 1990; Matthaeus and Zhou, 1989; Zhou and Matthaeus, 1990] and 
Dobrowolny et al.'s generalized Kraichnan phenomenology [Dobrowolny et 
al, 1980; Kraichnan, 1965]. The basic idea paper I is similar to that of Tu et 
al [1984] and Tu [1988] which is that the dissipation rates can be calculated 
from the inertial range energy spectra without any knowledge of the details 
of dissipation mechanism. There are some differences between our model and 
Tu's [1988] model; they are compared in paper 1. In paper 1 the constants of 
the phenomenologies were treated as free parameters. Their results were in 
general agreement with the observations for certain values of the free param- 
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eters. Recently Verma and Bhattacharjee [1995b, c] have calculated the con- 
stants of the Kolmogorov-like phenomenology theoretically using the direct 
interaction approximation (DIA) technique of Kraichnan [Kraichnan, 1959]. 
The constants calculated theoretically are in general agreement with the re- 
sults of numerical simulations by Verma [1994] and Verma et al. [1995d]. 
However, the constants calculated by DIA and numerical simulations arc dif- 
ferent from the ones used in paper I. In this paper we calculate the dissipation 
rates using the constants calculated theoretically, contrary to the procedure 
of paper I in which the constants were basically free parameters. Section 
2 of this paper contains these calculations. Note that the dissipation rates 
in the solar wind varies with distances. In this paper we estimate typical 
dissipation rates in the solar wind (at 1 AU). 

In section 3 we estimate the contributions of turbulent thermal conduc- 
tion to the heating in the solar wind. In section 4 of this paper we have 
estimated turbulent dissipation rates in the solar corona using Dobrowolny 
et al.'s generalized Kraichnan phenomenology. Here we also restate Holl- 
weg^s [1986] arguments in terms of the Kolmogorov-like MHD turbulence 
phenomenology. Section 5 contains the concluding remarks. 



2 Turbulent Heating in the Solar Wind 

In paper I we have derived a formula for the temperature evolution of the 
solar wind protons due to turbulent heating. Other source of heating were not 
included in that paper. However, one can easily generalize the equation for 
the temperature evolution when heating due to shocks and heat conduction 
are also included [see Priest, 1982]. The generalized equation is 

dT 2RT _ E ( V-q\ 

1" 7^ — ^turb + ^shock H ) (Ij 

dr Cy r UCv \ p J 

where T is the temperature, r is the radial distance from the sun, R is the 
Rydberg's constant, Cy is the specific heat per unit mass at constant volume, 
E is the mass per unit mole (~1 gm for the solar wind plasma), U is the solar 
wind mean speed, eturb is the turbulent dissipation rate per unit mass, eshock 
is the dissipation rate per unit mass due to shocks, q is the heat flux, and 
p is the density. Here as well as in paper I we have assumed constant speed 
for the solar wind, hence p cx r~^. We make this assumption to simplify the 
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analysis even though the density is observed to vary somewhat differently 
from r~^. The Eq. (P has the immediate solution 



-4/3 
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To{U) + 777^ /(f) ( eturb + ^ shock 




(2) 

where To(t/) is the temperature at ro for a stream with velocity U. 

In paper I we have calculated the turbulent heating in the solar wind using 
the Kolmogorov-like MHD turbulence phenomenology and Dobrowolny et 
al.'s modified Kraichnan's phenomenology. Here we state the energy spectra 
in these phenomenologies. In Kolmogorov-like phenomenology, they are 

E^{k) = C± (e±)'^' [e^y'^'k-'/^ (3) 

where E^{k) are the energy spectra of fluctuations = u±b (u is the 
velocity fluctuation, and b is the magnetic fluctuation in velocity units), 
are the dissipation rates of z^, and are Kolmogorov's constants for 
MHD turbulence. By inverting the above equation we can obtain the total 
dissipation rate which is 

E+ik)f'k'>/^ (4) 




where ai = E^/E^. Note that using this scheme we can obtain the dissipa- 
tion rates using the energy spectra of the wind and the constants without 
any knowledge of the details of dissipation mechanism. 

In this paper, without any loss of generality we take ai < 1. The normal- 
ized cross helicity dc, a measure of velocity and magnetic field correlation, is 
defined as 

2v ■ b _ 1 - ai 

The normalized cross helicity plays an important role in determining dissi- 
pation in the solar wind [Verma et ai, 1995a]. 

Recently Verma and Bhattacharjee [1995b, c] have calculated the values 
of the constants using DIA. They find that the are not universal 
constants as is the case in fiuid turbulence, but they depend on the Alfven 
ratio ta (ratio of kinetic and magnetic energy) and the normalized cross 
helicity o"c. In Table 1 we list some of the values of these constants. In these 



4 



calculations they take -E^ > E^. Unfortunately, in paper I it was assumed 
that = C~ = C, an approximation valid only for nonAlfvenic streams. 
In this paper we will use constants and obtained by DIA. 

In Dobrowolny et al.'s generalized Kraichnan's phenomenology, the dissi- 
pation rates are given by 

€+ = €- = A-^B^^E+{k)E~{k)k' (5) 

where A is Kraichnan's constant, and Bq is the mean magnetic field or the 
magnetic field of the largest eddies. Verma and Bhattacharjee [1995c] have 
also calculated the Kraichnan's constant for a variation of Dobrowolny et 
a/.'s model and found that the Kraichnan's constant A is approximately 

2.0 for nonAlfvenic streams. In this paper we emphasize the Kolmogorov- 
like phenomenology because its predictions appear most consistent with the 
observed energy spectra of —5/3 for both in the solar wind [Matthaeus 
and Goldstein, 1982; Marsch and Tu, 1989]. 

From the energy spectra and the constants we can estimate the dissi- 
pation rate Sturb at reference position tq using Eq. @). Then we can obtain 
the temperature at any position r using Eq. (^. In paper I the best fit to the 
observed temperature evolution of the solar wind was found when C = 1.0 
for nonAlfvenic wind and C = 8.0 for Alfvenic wind, and the dissipation 
rates at 1 AU was found to be of the order of 10"^ km^/s^. The calculation 
of the temperature evolution using Dobrowolny et a/.'s model with neither 
A = 1 nor A = 0.5 provided a good fit to the observed temperature evolu- 
tion, however here too for A = 1, the dissipation rates at 1 AU were of the 
order of 0.5 x 10~^ km^/s^ [Verma et ai, 1995a]. See paper I for details. 

The turbulent dissipation rate depends on the constant C, and it is crucial 
to use the correct constants C for proper estimation of turbulent heating. 
In the following subsections we use the constants calculated by Verma and 
Bhattacharjee [1995b, c] to estimate the turbulent dissipation rates in the 
solar wind. 

2.1 NonAlfvenic Streams 

In the outer heliosphere most of the solar wind streams have been observed 
to be nonAlfvenic (cTc ~ 0) with ~ 0.5. From Table 1 C = 4.04 [Verma 
and Bhattacharjee, 1995b, c] for these streams. In paper I Verma et al. found 
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that when C = 1.0, turbulent heating can account for all the heating in 
the solar wind. They estimated turbulent dissipation rate at 1 AU when 
C = 1 to be approximately 10~^ km^s~^. However, since theoretical value 
of C is approximately 4 in the outer heliosphere [Verma and Bhattacharjee, 
1995b, c], the correct turbulent dissipation rate is approximately (refer to Eq. 

(i) 

6 = 6„/ (^il^j - 1-25 X \^-'km^s-\ (6) 

Hence, in the outer heliosphere, turbulent heating is not sufficient to heat 
the plasma to the observed temperature. The rest of the heating can be 
provided by shocks, thermal diffusion, and other mechanisms. Whang et 
al. [1990] showed that the shock heating is significant, while Gazis [1984] 
showed that heating due to thermal conduction is significant (for thermal 
conduction, also refer to Marsch et al. [1983]; Marsch and Richter [1984]). 
We briefiy sketch the results of their work in the section 3. 

In the inner heliosphere, the fiuctuations are somewhat fiuid dominated. 
Therefore, C in the inner heliosphere is smaller than C of the outer helio- 
sphere (see Table 1), hence resulting in a larger turbulent heating. If we take 
ta in the range of 1 to 2, C is in the range of 3.38 to 2.58 (see Table 1). 
With these constants too the turbulent heating is utmost C~^^'^ ~ 25 — 30% 
(approximately one-quarter) of the of observed heating (heating rate of pa- 
per I) in the solar wind. The corotating shocks are absent in the inner 
heliosphere. Therefore, shock heating is presumably not significant for non- 
Alfvenic streams in the inner heliosphere. Hence, neither shocks nor turbu- 
lence can provide enough heating in the inner heliosphere. In section 3 we 
will estimate the contributions of thermal conduction to the overall heating 
in the inner heliosphere. 

Regarding Dobrowolny et a/.'s model, for nonAlfvenic streams, the Kraich- 
nan's constant A is approximately 2.0. With A = 2.0, the dissipation rates 
at 1 AU will be of the order of epaperi/A"^ ~ 0.5 x lO^V^^ ~ 10"^ km^s'^ 
(see Eq. [Q), which again is lower than the observed heating. 

2.2 Alfvenic streams 

In paper I Verma et al. find that the observed temperature evolution of the 
Alfvenic streams are in good agreement with the predictions of Kolmogorov- 
like phenomenology if Kolmogorov's constant for MHD turbulence C is equal 
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to 8.0. In paper I and were taken to be equal which is inconsistent 
with the DIA calculations of Verma and Bhattacharjee [1995b, c]. In the 
following discussion we modify the dissipation rates of paper I using the 
theoretical results of Verma and Bhattacharjee [1995b, c]. 
For large cTc, since ai is small, from Eq. (H) 



If we take dc = 0.90 (ai = 0.05) and ta = 2.0, Verma and Bhattachar- 
jee [1995c] predict = 1.54 and = 2.08. With these constants eturb 
predicted by the Kolmogorov-like phenomenology is larger by a factor of 
g3/2^ (C^V C^) ~ 10 as compared to what was calculated in paper I, i.e., 
the predicted tturb with C"^ of DIA is 10~^ km^s~^. When we take = 1.0, 
unfortunately, the maximum cTc we can get in DIA scheme is 0.5 (after Uc 
= 0.5 there is no consistent solution). This limitation indicates that some of 
the assumptions of Verma and Bhattacharjee [1995b, c] are not fully consis- 
tent. However, from the general trend of we believe that the Kolmogorov- 
like MHD turbulence phenomenology predicts higher dissipation rate as com- 
pared to what is observed in the solar wind. 

The above study appears to justify the conjecture of Marsch [1991] and 
Grappin et al, [1991] which states that in the solar wind the Alfvenic streams 
have not reached steady-state, and the turbulence in these streams is not 
yet fully developed. The nonlinear interactions may be distributing energy 
among large and intermediate wavenumbers, with only a small amount of 
energy flowing into the dissipation range. Therefore, the real dissipation rate 
is lower than that predicted by the phenomenology, which assumes that the 
system is in steady-state. Any modelling of "approach towards steady-state 
turbulence" will shed light to the evolution of energy spectra and help us in 
correct estimation of turbulent heating in the Alfvenic streams in the solar 
wind. 

The simulational studies show that the system reaches steady-state after 
around 1 eddy turnover time [Verma, 1994; Verma et ai, 1995d]. For solar 
wind it corresponds to the distance of approximately 2 AU {l{l/urms)Vsw = 
1 * {O.lAU/20) * 400 = 2AU). Hence, it is possible that the Alfvenic streams 
have not reached steady-state during its transit in the inner heliosphere, 
and the estimated turbulent heating by the turbulence phenomenologies is 
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higher than reahstic heating in the solar wind. However, it is hard to jus- 
tify the above mentioned conjecture in the outer hehosphere. Also note the 
nonAlfvenic streams appear to have reached steady-state by the time they 
reach 1 AU. This could be due to the initial conditions of the nonAlfvenic 
streams near the sun. We need to perform careful analysis of the approach 
to steady-state and effects of cross helicity in MHD turbulence to understand 
this behaviour. 

Considering that the turbulent heating is not sufficient to heat the non- 
Alfvenic streams to the observed temperature, in the following section we 
consider the heating due to shocks and thermal diffusion. 

3 Heating due to Shocks and Thermal Diffu- 
sion 

Whang et al. [1990] studied the average entropy increase across a shock 
for heliocentric distances of 1-15 AU. They showed that an average entropy 
increase per shock, AS", is approximately 0.8 x 10"^^ J/K/proton. The 
strong shocks were in the region of 1-5 AU. Whang et al. 's [1990] found 
approximately 400 shocks in 1-15 AU region. Therefore, the dissipation rate 
^shock = NT AS/ At, where N is the number of shocks in 1-15 AU, T is the 
temperature and At is the time taken to travel 1 AU to 15 AU, will be 

eshock NT AS/ At 

~ 400 X 10^ X 0.8 X 10-23/(14 x 1.5 x 107500)J/proton 

~ 400 x 10^ x 0.8 X 10-23 X 10-V(mp x 14 x 1.5 x 10V500)km7s^ 

~ lO-^kmVs^ 

(8) 

Here nip is the mass of a proton. From the above calculation we see that 
shock heating is comparable to the turbulent heating, and they are one of the 
major sources of heating in the outer heliosphere. However, note that Whang 
et al. 's [1990] observational study of entropy increase in the solar wind with 
the radial distance 

S-Se^ 3.81 X log(r/rE) 10-23 J/K/proton (9) 

where the subscript E denotes the conditions at rg = 1 AU, is the total 
entropy increase due to all the heating sources present in the solar wind, not 
just due to shocks, as inferred in the Whang et a/.'s [1990] paper. 
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Regarding heat conduction, it is widely believed that the coupling be- 
tween electrons and protons is weak because the protons-electron collision 
frequency is much smaller as compared to electron-electron or proton-proton 
collision frequencies. Due to this reason, for the following discussion we as- 
sume that the temperature of the protons in the solar wind are affected by 
the heat flux of protons alone, not by the heat flux of electrons. Marsch and 
Richter [1984] performed observational studies of the solar wind at 0.95 AU 
and reported that the heat conduction by protons qi is speed dependent and 
is approximately 10~^ ergs/cm^s (also refer to Hundhausen [1972], Marsch 
et al. [1983], and Gazis [1984]). At 0.95 AU, using particle density as 5/cc, 
the proton heating rate Vqi/ p due to this heat flux is 10~^ km^s~^, which is 
comparable to the heating rates due to turbulence. Hence, the dissipation 
due to heat flux appears to be a likely candidate for the solar wind heating. 

Now we will obtain theoretical estimates of the heating due to thermal 
conduction using order of magnitude calculations and compare them with 
the observational results. The source of heating due to thermal conduction, 
according to Eq. ([^), is 

i/ = = ^i^^ . £lZi, (10) 

p p r^p 

where Ki are the thermal conductivity of protons, are the temperature of 
protons, and r is the radial distance from the sun. Here we take r = 1 AU 
and Tj = l^f'K. Note that the quantities H decreases with distance. The 
proton thermal conductivity ~ pcyvl ~ pCyV [Priest, 1982], where cy is 
the specific heat per unit mass of protons, v is the relevant velocity scale, I 
is the length-scale, and v is the kinematic viscosity. Therefore, 

(11) 

Note that R = 8.31 J/K/mole, cy for protons is 12.5 Joule/K/gm. 

Montgomery [1983] estimated viscosity in the solar wind using the clas- 
sical transport theory of Braginskii [1965] and obtained typical u ~ 10^^ 
km^s~^ [Montgomery, 1983]. We assume that the proton temperature are 
approximately 10^ K. With the above viscosity we find that at 1 AU (r = 1 
AU), H ~ 10~^^ km^s~^. This H is negligible as compared to the turbu- 
lent dissipation rates and is also inconsistent with the H obtained from the 
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observations [Marsch and Richter, 1984; Marsch et ai, 1983; Hundhausen, 
1972]. 

However, since the solar wind plasma is turbulent, thermal conduction 
coefficient should be determined by large-scale velocity and length scales 
[Landau, 1987]. The basic idea is that in a turbulent plasma, the large scale 
eddies carry heat flux from one region to the other. For this reason, we must 
use turbulent eddy viscosity due to the large eddies in our calculation for 
K. Therefore, v ~ vl, where f ~ 20 km/s (speed of large eddies) and I 
~ 10^^ cm (size of the large eddies). This yields H ~ 10~^ km^s~^. Since our 
calculations involve only order of magnitude estimates, we can say that our 
estimates are in general agreement with the observational results of Marsch 
and Richter [1984], Marsch et al. [1983], and Hundhausen [1972]. Therefore, 
heating due to thermal conduction, corotating shocks, as well as turbulence 
appear to be the heating sources of the solar wind in the outer heliosphere. 

Gazis [1984] estimated the thermal energy flux at 1 AU in the solar wind 
and found it to be equal to (2.5 ± 1.0) x 10~^ ergs/cm^s, which corresponds 
to V ■ q/p ~ 2 X 10~^ km^s~^, that is approximately 10 times larger than the 
dissipation rate estimated in paper I. However, note that in Gazzs' formalism 
thermal condition was considered as the only source of entropy increase (see 
Eq. (24) of Gazis [1984]), therefore, the heating rate calculated by Gazis 
was an over-estimate. This is not surprising since Gazis [1984] sets an upper 
limit on the total heat flux q. 

Regarding heating due to thermal conduction in the inner heliosphere, 
the substitution of = lO^iT, r = 0.3 AU, v = km/s, and / = 10^^ 
Eq. ([TT|) yields Hi ~ 10~^ km^s"'^. This quantity is comparable to turbulent 
heating rate in the solar wind. Thus, ion heating due to thermal conduction 
is a potentially significant heating source in the inner heliosphere and could 
provide significant part of the heating in the inner heliosphere. Note that the 
corotating shocks are absent in the inner heliosphere, therefore, shock heating 
is presumably not significant for nonAlfvenic streams in the inner heliosphere. 
Of course, other sources, e.g., stream-stream interactions, neutral ions etc. 
may also provide significant heating in the inner heliosphere. 

To summarize, turbulence provides only part of the dissipation in the solar 
wind. Corotating shocks as well as thermal conduction appear to be other 
heating sources in outer heliosphere and could provide rest of the heating. 
In the inner heliosphere, heating due to thermal conduction appear to be 
significant. In absence of corotating shocks in the inner heliosphere, thermal 



10 



condTiction may provide significant part of tlie dissipation here. In this paper 
we have performed only order of magnitude estimates of these quantities, 
therefore, we cannot make any definite statements regarding energy budget 
in the solar wind. 

Turbulent heating is considered to be one of the potential mechanism 
which could explain solar corona heating problem. In the following section 
we estimate amount of turbulent heating in the solar corona using the above 
mentioned MHD turbulence phenomenologies. 

4 Turbulent Heating in the Solar Corona 

The temperature of the sun decreases radially outwards to about 6000 K 
at photosphere, then it rises abruptly upto 2 — 3 x 10^ K at the corona. 
Why the coronal temperature rises so steeply remains an open problem. 
Various heating mechanisms have been proposed, turbulent heating being 
one of them [Hollweg, 1984; Browning and Priest, 1984; Browning, 1991 
and references therein]. In the following discussion we state Hollweg^s [1984] 
calculation of the dissipation rates in the solar corona using the Kolmogorov- 
like phenomenology, and we calculate the dissipation rates using Dobrowolny 
et al.'s generalized phenomenology. We compare these estimates with the 
observational results. Here too we use dimensional arguments similar to 
those used for solar wind. The power of these arguments lies in the fact 
that we can estimate the dissipation rates from the inertial range spectra 
or from the energy fed in at large-scales, but without any knowledge of the 
dissipation mechanism which may be active at small-scales. 

Here we rephrase the Hollweg 's arguments [1984] on the estimation of 
turbulent heating in the corona in terms of E^{k). We assume that in the 
corona E^{k) ~ E~(k) ~ E{k). Therefore, the dissipation rate according to 
the Kolmogorov-like phenomenology will be given by 

e ~ C-^/''{E{k)kf'''k. (12) 

Taking the average speed of the large-scale eddies as 30 km/s, E{k)k ~ 
v1 ~ 10^ km^s~^. The size of the largest eddies is roughly 10^ km. Therefore, 
e ~ 10 krn^s"^. Hollweg [1984] argues that the above dissipation rate is 
compatible with the heating requirements of coronal active region loops. 
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In contrast, the dissipation rate in Dobrowolny et al.'s generalized phe- 
nomenology [Dobrowolny et al, 1980] will be 

e^A-^B^\E{k)kfk. (13) 

With mean number density as 10^^ particles/m^ and the magnetic field of the 
order 10-100 mT, Bq in velocity units is 10^ — lO'' km/s. Taking E{k)k ~ 10^ 
km^s^^ as before, we obtain e ~ 10^^ — 10^'^ km^ s^^. This dissipation 
rate is three to four order of magnitude lower than that estimated by the 
Kolmogorov-like model. Therefore, if turbulent heating is playing a major 
role in the coronal heating, the Kolmogorov-like phenomenology must be at 
work in the solar corona. This is surprising because one would expect Do- 
browolny et al.'s model to work due to the dominance of the mean magnetic 
field over the fluctuations. In the solar wind too, even though the Parker 
field is stronger than the fiuctuations, the energy spectra follows roughly 
k~^/^ power law. Also the temperature evolution of the solar wind is better 
explained by the Kolmogorov-like model rather than Dobrowolny et al.'s phe- 
nomenology (see paper I). These observations indicate that Kolmogorov-like 
phenomenology is applicable under conditions where it is not expected to 
work. 

The eddy turnover time-scale in the solar corona is l/vi ~ 10^/30 ~ 300 
sec. Since the plasma reaches steady-state in around 1 eddy-turnover time 
[Verma, 1994; Verma et al., 1995d], the question arises whether it is possible 
for the plasma to reach steady-state and be able to supply enough energy 
in the dissipation range. It is suggested that on the solar surface and in 
the coronal loops plasma could stay for around 300 seconds, and therefore, 
turbulent heating could be one of the possible candidate for coronal heating 
(see Gomez and Fontdn [1988]; Browning [1991] and references therein). To 
resolve the above issues one needs to resort to model calculations; for some 
of the recent calculations refer to Gomez and Fontdn [1988], Heyvaert and 
Priest [1992], ZirA;er [1993], Vekstein et al. [1993]. 

5 Conclusions 

In this paper we have estimated the contributions of turbulent dissipation 
in the heating of the solar wind using Kolmogorov's constants for MHD tur- 
bulence calculated theoretically by Verma and Bhattacharjee [1995b, c]. We 
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find that for tlie nonAlfvcnic streams, turbulence contributes only partly 
to the total heating in the solar wind. The remaining part of the heating 
should be provided by corotating shocks, thermal conduction, stream-stream 
interactions, interactions with the neutral ions, and by other sources. The 
corotating shocks are present in the outer heliosphere. As shown by Whang 
et al. [1990], they could be one of the major heating sources in the outer 
heliospliere, at least in 1-15 AU where they have been studied, and could 
provide a major fraction of the heating. The order of magnitude calculations 
show that the heating due to thermal conduction is significant both in the 
inner as well as outer heliosphere. In absence of corotating shocks in the 
inner heliosphere, thermal conduction probably plays an important role in 
the solar wind heating and could possible provide significant fraction of the 
heating here. Our results are based on order of magnitude estimates, there- 
fore, we cannot make definite statements. Models incorporating these results 
may help us in making definite predictions about the contributions by these 
sources to the solar wind heating. 

For the Alfvenic streams, the Kolmogorov-like MHD turbulence phe- 
nomenology with Kolmogorov's constants calculated by Verma and Bhat- 
tacharjee [1995b, c] predicts higher dissipation rate as compared to what 
is observed in the wind. The resolution of the paradox could be that the 
Alfvenic streams have not reached steady-state [Marsch, 1991 and references 
therein; Grappin et al., 1991], and the energy is just being distributed among 
various modes, with only a small amount of energy flowing into the dissipa- 
tion range and heating the plasma. Hence the dissipation rates predicted 
by the Kolmogorov-like MHD turbulence phenomenology are over-estimates. 
The studies on "approach to steady-state" will help us in proper estimation 
of turbulent heating in the Alfvenic streams and in understanding of other 
related problems. 

Lastly, we have discussed the contributions of turbulent dissipation in 

the coronal heating. Hollweg [1984] has calculated the heating rate using 
Kolmogorov's fiuid turbulence phenomenology and showed that the results 
are consistent with the required energy flux in the coronal active loops. In 
this paper we calculate the heating rates using Dobrowolny et al. 's general- 
ized model and show that the prediction from this model is four orders of 
magnitude smaller than the required heating rates. This calculation shows 
that if turbulent heating is major source of heating in the corona, then 
the Kolmogorov-like phenomenology rather than Dobrowolny et a/.'s phe- 
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nomenology is valid in the corona. In the corona, the mean magnetic field 
in velocity units is much larger than the fiuctuations, hence, if we assume 
that the heating in the solar corona is due to turbulence, the this result is 
contrary to the assumptions of the phenomenologies. However, the result is 
in agreement with the solar wind results, where the Kolmogorov-like model 
works even though the Parker's field is larger than the fiuctuation. This is an 
important puzzle in MHD turbulence, and future theoretical and numerical 
work will help in clarifying this issue. 

The author thanks D. A. Roberts, M. L. Goldstein, and J. K. Bhattachar- 
jee for discussion and comments. 
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rA 




C+ 


C~ 


0.5 


0.0 


4.04 


4.04 


1.0 


0.0 


3.38 


3.38 


2.0 


0.0 


2.58 


2.58 


5.0 


0.0 


1.92 


1.92 


1.0 


0.43 


0.54 


53.80 


1.5 


0.67 


1.44 


6.47 


2.0 


0.90 


1.54 


2.08 



Table 1: Kolmogorov's constants for MHD turbulence for various values 
of Alfven ratio and normalized cross helicity [Verma and Bhattacharjee, 
1995b,c] 
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